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m ■ ABSTRACT 

Among the most remarkable features of the stellar population of R136, the central, 
young, massive star cluster in the 30 Doradus complex of the Large Magellanic Cloud, 
are the single stars whose masses substantially exceed the canonical stellar upper 
mass limit of 150M Q . A recent study by us, viz., that of Banerjee, Kroupa & Oh 
. (2012; hereafter Paper I) indicates that such "super-canonical" (hereafter SC) stars 

Q-i can be formed out of a dense stellar population with a canonical initial mass func- 

L , tion (IMF) through dynamically induced mergers of the most massive binaries. The 

above study consists of realistic N-body computations of fully mass-segregated star 
' clusters mimicking R136 in which all the massive stars are in primordial binaries. In 

C^j . the present work, we study the formation of SC stars in the computed R136 models 

of Paper I in detail. Taking into consideration that extraneous SC stars form in the 
computed models of Paper I due to the primordial binaries' initial eccentricities, we 
compute additional models where all primordial binaries are initially circular. We also 
• take into account the evolution of the mass of the SC stars and the resulting lifetime 

^vq , in their SC phase using detailed stellar evolutionary models over the SC mass range 

00 ■ that incorporate updated treatments of the stellar winds. In all these computations, 

1 we find that SC stars begin to form via dynamical mergers of massive binaries from 

. w 1 Myr cluster age. We obtain SC stars with initial masses up to « 250Af© from 

j^-*) ■ these computations. Multiple SC stars are found to remain bound to the cluster si- 

| multaneously within a SC-lifetime. However, we also note that SC stars can be formed 

■ a t runaway velocities which escape the cluster at birth. These properties of the dy- 

r* ■ namically formed SC stars, as obtained from our computations, are consistent with 

, \ the observed SC stellar population in R136. In fact, the evolutionary models of SC 

. stars imply that had they formed primordially along with the rest of the R136 cluster, 

5^ 1 i.e., violating the canonical upper limit, they would have evolved below the canonical 

150A/ Q limit by w 3 Myr, the likely age of R136. Thus according to the new stellar 
evolutionary models, primordially-formed SC stars should not be observable at the 
present time in R136. This strongly supports the dynamical formation scenario of the 
observed SC stars in R136. 

Key words: methods: numerical - stars: kinematics and dynamics - stars: luminosity 
function, mass function - stars: massive - stars: winds, outflows - open clusters and 
associations: individual (R136) 



1 INTRODUCTION 

The study of the functional form of the number distri- 
bution of stars in galaxies, with which they are born, or 
the stellar initial mass f unction (IMF) has always been of 
fundamental importance (|Bastian et al.ll2~010l ; iKroupa et al.l 



* E-mail: sambaran@astro.uni-bonn.de (SB) 
f E-mail: pavel@astro.uni-bonn.de (PK) 
J E-mail: skoh@astro.uni-bonn.de (SKOH) 



l2012t ). The high-mass end of the stellar IMF is in partic- 
ular focus due to the feedback it gives to the star-forming 
gas in the form of radiation pressure, kinetic energy from 
the stellar wind and supernova ejecta and as well the chem- 
ical enrichment from the latter source. An upper limit of 
~ 6OM0 to the mass of a star could be set by the "Ed- 
dington limit" l|Eddingtonl Il926t) or the point of balance 
of gravity by the radiation pressure of a star. However, 
in practice stellar masses can easily exceed the Edding- 
ton limit since massive stars are not fully radiative but 
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contain convective cores (|Kippenhan fc WeigertH"l990l ). The 
next upper limit is set by the possi bility of the destruction 
of a s t ar due to thermal puls at ions llSchwarzschild fc Harm! 
1 19591 ; iBeech fc Mitalasl ll994l V IStothers! l| 19921 ) determined 
this limit to be at m max , « 120 - 15OAf for [Fe/H]« and 
m max „ « 9OM for [Fe/H]« -1. 

Difficulties also arise if one considers the growth of a 
proto-star via gas accretion which is crucial for massive star 
formation. Here one again encounters the w 6OAf Edding- 
ton limit. Stellar formation models lead to a mass limit near 
40— 1OO M imposed by feedback on a spherica l accretion en- 
velope (|Kahnlll97llWolfire fc Cassinellill 19871 ) . Some obser- 
vations suggest that proto-stars may be accreting material 
in disks rather than spheres (e.o.. lChini et al]|2004l ) in which 
case it may be possible to overcome the radiation pressure 
at the equator of the proto-star. Studies on the formation 
of massive stars through disk-accretion with high accretion 
rates, thereby allowi ng the radiation to esc ape preferentially 
along the poles (e.a. Jjiuna fc Adams! 199r3 ) indeed allow for- 
mation of stars with larger masses. 

The feedback-induced mass limit can al so be avoided 
if massive stars can form thro ugh mergers dBonnell et al.l 
1998|; IZinnecker fc Yorkei |2007| ). In this scenario, massive 
stars form via coalescence of intermediate-mass proto-stars 
in the cores of dense stellar clusters that have undergone 
core-contraction due to rapid accretion of gas with low spe- 
cific angular momentum. The high central density required 
for the mergers (10 s Af pc -3 ) is still difficult to achieve but 
it should be noted that an observable young cluster is nec- 
essarily disposed of a substantial fraction of its natal cloud 
and is thus likely to be always observed in an expanding and 
hence diluted phase. Recently-set limits on the radii of em- 
bedded clusters indeed suggest them to form very compact 
l|Marks fc Kroupall2012T l. In this context, it is worthwhile to 
note that in the computations discussed in the following sec- 
tions we do get binary coalescence events even in moderately 
dense stellar clusters. This happens due to the presence of 
a substantial number of binaries in the cluster core which 
have much larger encounter cross sections than single stars 
due to their bigger geometrical sizes. 

While the existence of an upper limit of the stellar 
mass is rather obvious from theory since several decades, 
such an upper limit has been established from observations 
only recently. There have been some earlier indications of 
the presence of a cut-off near m max * ~ 150Mq from ob- 
servations of young massive clusters, particularly in R136 
|Massev fc Hunted Il998l; lMassevll2003h in the LMC and in 
Arches cluster (|Figerll2003l ) close to our Galactic center, but 
these results were not sufficiently convincing in the sense 
that no statistical significance were attached to them. The 
observed upper limit was consider ed to be a limitation due to 
sampl ing rather than a true limit l|MassevH2003h . lElmegreenl 
(2000) also noted that random sampling from an unlimited 
IMF for all star-forming regions in the Milky Way would 
lead to the prediction of stars with masses > IOOOA/q, un- 
less there is a sharp down-turn in the IMF beyond several 
1OQM . 

IWeidner fc Kroupal l|2004h for the first time gave a criti- 
cal look at the question of the existence of a physical upper- 
limit of w 15OA/ in the IMF, as observed in R136 (the 
central massive star cluster in the 30 Doradus region of the 
LMC). Assuming R136 has a mass in the range 5 x 10 4 < 



M c i < 2.5 x 1O 5 M (jSehnan et al.lll999l L they found that for 
a canonical IMF with m max * = oo, 10 < N(> 15OM ) < 40 
stars are missing in R136 that have masses > 15OM . The 
probability that no stars are observed among the 10 ex- 
pected ones, assuming m max * = oo, is p = 4.5 x 10~ 5 , i.e., 
the observed massive stellar content of R136 implies a phys- 
ical st ellar mass limit at m max » w 15OM . Similarly, iFigerl 
(|2005h found a dearth of N(> 15OAf ) = 33 stars in the 
Arches cluster, where the shallowing of its stellar mass func- 
tion due to its rapid tidal dissolution has been incorporated. 
The corresponding number of missing stars for a canonical 
IMF is JV(> 15OM ) = 18 which gives a non-detection prob- 
ability as small as p = 10 -8 . Thus the Arches cluster also 
has m max « 1 5OAf with a very high significance. Finally, 
lOev fc Clarke] (|2005l ) studied 9 clusters and associations in 
the Milky Way, LMC and SMC to investigate the expected 
masses of the most massive stars in these for different up- 
per mass limits (120, 150, 200, 1000 and 10000M©). They 
concluded that the observed number of massive stars sup- 
ports the existence of a general physical upper mass cutoff 
within the range 12OM < m max * < 2OOM with a high 
significance. 

The massive stellar population of young massive star 
clusters therefore indicates a physical upper mass limit near 
?rim ax * ~ 15OAf which is the canonical upper limit of the 
stellar IMF. One can regard this mass as the limit imposed 
by the process of star formation under physical conditions 
achievable in a star-forming core. The origin of this very 
limit is under investigation and the v alue should currently be 
taken as empirical only. As coined bv lKroupa et al.l l|2012l ). a 
stellar population containing stars with their zero-age- main- 
sequence (ZAMS) masses up to 15OM is "saturated". 

In this paper, our concern is again related to the massive 
stellar population in R136. The part icular aspect tha t we fo - 
cus on is related to a recent study bv lCrowther et ail (|2010T ). 
These authors re-analysed the massive stellar population of 
R136 in unprecedented detail using Hubble Space Telescope 
and Very Large Telescope spectroscopy and high spatial res- 
olution near-IR photometry to find 4 stars, within the cen- 
tral 1 x 1 pc of R136, with masses 165 — 32OM , i.e., sub- 
stantially above the canonical limit One can call such a stel- 
lar population, containing single stars with masses substan- 
tially exceeding 15OM , as "super-saturated" and the stars 
surpassing the canonical upper limit can be called "super- 
canonical" (|Kroupa et al.|[2012l ; hereafter SC). Although the 
stellar population of R13 6 has been studi e d by earlier au- 
thors, the crucial turn in ICrowther et al.l |2010) is due to 
these authors' consideration that the observed SC stars in 
R136, in spite of exhibiting WN-type spectra and possess- 
ing strong winds, are actually young, main-sequence stars 
rather than cl assical Wolf-Rayet (W R) stars. To take this 
into account, ICrowther et alj (|2010l ) incorporate detailed 
stellar evolutionary models of very massive main-sequence 
stars that include state-of-the-art treatment of stellar winds 
which lead to the inference of the SC masses. 

More recently however, iBaneriee. Kroupa fc Ohl (|2012T ) 
(hereafter Paper I) showed that single stars with masses 
m s > 15OM can indeed form in R136 through mergers 
of the members in massive binaries. These authors com- 
puted the evolution of model clusters that resemble R136 
using the st ate-of-the-art direct N-body integration code 
"NBODY6" (|Aarsetbll2003h . In these computationally chal- 
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Figure 1. Formation of SC stars in the computed "Model 3" (upper panel) and "Model 4" (lower panel) of iBaneriee. Kroupa fc Ohl 
(2012). The panels show the masses of the SC stars at the times of their first appearance and their subsequent mass depletion due to 
stellar wind as obtained from within NBODY6. For "Model 3" (upper panel), the SC stars represented by the grey lines, that appear 
from the very beginning, are the "spurious" ones in this example in the sense explained in Sec. 12.31 Similar description applies for the 
lower panel. 



lenging models, all the O-stars w ere taken to be in tight 
binaries as observations indicate l|Sana fc Evans! l201ll ; SCO 
Sec. [2j but the ZAMS mass of a single component never ex- 
ceeds m max * = 15OA/0. They found that in course of their 
dynamical evolution, the massive binaries in the model clus- 
ters can often merge due to hardening and/or eccentricity 
enhancement due to dynamical encounters to produce single 
stars of m s > 15OM0. This implies that the presence of SC 
stars in R136 does not illustrate a violation of the canonical 
upper limit. 

Of course, the value of the u pper stellar mass lim it is not 
necessarily exactly = 15OM0 (|Oev fe Clarke! 120051 ). How- 
ever, for definiteness, we take the following hypothesis: there 
exists a fundamental upper mass limit of m max * = 15OM0 or 
the canonical limit of ZAMS stars formed in clusters and any 
observed more massive ("super-canonical") star is created 
from dynamically-induced stellar mergers within the dense 
young cluster. In this study, we aim to test the above hy- 
pothesis or answer the following question: Can the observed 
number of super-canonical stars in R136 be explained by the 
naturally occurring stellar- dynamical and stellar- evolution 
processes in R136? To that end, we utilize the 4 direct N- 
body computed massive cluster models of Paper I mimicking 
R136 (Sec. l2.ll fc l2.2p and additionally perform 5 more simi- 
lar computations (Sec. 12. 3) ) to trace the formation of SC stars 
through dynamical means. These computed models initiate 



with properties that are consistent with the observed prop- 
erties o f young clusters, in particular, primordi al mass segre- 
gation jLittlefair et al.ll2003l ; IChen et alj|2007h and massive 
stellar binaries with component mass ratio close to unity 
ijSana fc Evansl [201 lh . Furthermore, we utilize detailed ro- 
tating stellar evolution models in the SC mass range to esti- 
mate the mass evolution of the dynamically formed SC stars 
and the resulting lifetimes in their SC phases (Sec. Q. We 
conclude the paper by summarizing our results and pointing 
out the limitations of the present study (Sec. [5J. 



2 COMPUTATIONS 

State-of-the-art calculations have been performed in Paper 
I to evolve model star clusters that mimic R136 in terms of 
observable parameters. The primary objective of these cal- 
culations was to study the ejected OB-stars from R136, in 
particular, whether massive runaway stars like VFTS 682 
and 30 Dor 016 ca n indeed be ejecte d dynamically from 
R136 as suspected l|Evans et al.l |2010| ; iBestenlehner et"al] 
120111 ). Besides good agreement of the kinematics of the mas- 
sive ejecta with those of these noted runaways, it is also 
noted in Paper I that the tight massive binaries which drive 
these runaways can merge due to the dynamical interactions 
leading to the formation of SC stars in each of the models, 
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Figure 2. Formation of SC stars in the models "C2", "C5" and "CIO" with all their primordial binaries initially in circular orbits (see 
Sec. I2.3( l. Unlike the computations in Paper I (c.f. Fig. [lj , there are no "spurious" SC stars in these cases. 



in agreement with the observed stellar population in R136 
|Crowther et alj|2010h . In this study, we continue to utilize 
the results of these computations. The initial conditions and 
the method of these calculations are accounted in detail in 
Sec. 2 of Paper I which we briefly recapitulate here. 

2.1 Initial conditions 

The above computations c omprise direc t N-body integra- 
tions of Plummer spheres l|Kroupall200ct ) with parameters 
conforming with those of R136. The initial mass of the Plum- 
mer sphe res is M e ;(0) m W 5 M (^ which is an upper limit 
for R136 (jCrowther et alJl201Ch and the half mass radii are 
taken to be rfc(O) « 0.8 pc. The clusters are made of stars 
with ZAMS masses drawn from a canonical IMF (|Kroupal 



l200lj ) over the range O.O8M0 < m s < 150A/q and are of 
metallicity Z = 0.5^0. 

As for the primordial binary population, stars with 
m s > 5Mq are all in binaries while all lighter stars are 
kept single. This termination of the binary population is 
due to computational ease; binaries bottleneck the calcu- 
lation speed of direct N-body integration significantly so 
that adopting a full spectrum of primordial binaries (i.e., 
100% initial binary fraction) in models of the size that we 
compute becomes prohibitive. Disregarding the binary pop- 
ulation for m s < 5Mq of course does not affect the ejec- 
tion of massive stars and mergers of massive binaries sig- 
nificantly (see Paper I for details). F ollowing the observe d 
period distribution of O-star binaries (jSana fc Evans|[201ll ). 
the orbital periods of the binaries having primary masses 
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m 3 > 2OM0 are chosen from a uniform distribution between 
0.5 < log 10 P < 4, where P is the orbital period in days. 
For primary m asses 5Mb < m s < 20M©, the "birth period 
distribution" of IKroupal (|l995l ) is ado pted. The eccentricity 
distri bution is chosen to be thermal l|Spitzerlll987l ; IKroupal 
2008). These dist ributions, how ever, are not modified via 
"eigenevolution" l|Kroupalll995h as the binary-systems' in- 
ternal interactions between massive pre-main-sequence stars 
are not yet quantified. 

The initial configu rations are fully mass-segregated 
(Bau mgardt et aH l2008) so that the massive binaries con- 
centrate within the clusters' central region. In other words, 
it is assumed that the clusters are primordially segre- 
gated w hich was the case for several Galactic globular 
clusters (iMarks fc Kroupal |2010T) a nd yo ung star clusters 
jLittlefair et all 120031 ; IChen et ail l2007t ). In these mass- 
segregated initial models, the stars are re-distributed over 
the cluster in such a way that the individual mass groups 
are in equipartition whilst the overall cluster is in virial equi- 
librium as well. The algorith m for creating such syste ms is 
described in the appendix of Baumga rdt et al.l (|2008T ). We 
emphasize that we are conforming with the observed states 
of massive, young clusters by starting our models binary-rich 
and mass-segregated. That is, we are not choosing some par- 
ticular fine-tuned initial conditions in order to enhance our 
results. 



2.2 N-body integration, stellar evolution and 
stellar mergers 

With such an initial setup, these computations are the very 
first direct N-body calculations under such extreme condi- 
tions (all massive stars in binaries, fully mass-segregated). 
To perform the direct N-body c omputations, the state- 
of-the-art "NBODY6" integrator (|Aarsethl 12003) has been 
used. In addition to integrating the particle orbits using 
the highly accurate fourth-order Hermite scheme and deal- 
ing with the diverging gravitational forces in close encoun- 
ters through regularizations, NBODY6 also employs the 
well-tested analytical stel lar and binary evolution recipes of 
iHurlev et ail (|2000l . 120021 ) (the SSE and the BSE schemes). 
These prescriptions are based on m odel stellar evolution 
tracks computed bv lPols et al.l l)l998T ). The wind mass loss of 
the massive main-s equence (hereafter MS) stars is ado pted 
using the empirical iNieuwenhuiizen fc de Jaeerl l|l990h for- 
mula. The mass loss rates are, of course, appropriately 
modified on the giant br anches and other e volved phases 
for stars of all masses ifHurlev et all l2000h . Admittedly, 
the INieuwenhuiizen fc de Jagen (jl99Clh mass loss rate is 



too simplistic for massive stars; it is likely to grossly un- 
derestimate the wind mass loss of massive stellar merger 
products, which we focus on here, as these objects usu- 
ally become Wo lf-Rayet (WR) stars even in their main 
sequence phase ijGlebbeek et al.l 120091 ). The above stellar 
wind recipe is currently the best available one in a direct 
N-body evolution code; other widely-used direct N-body 
codes, e.g., the "STARL AB" JPortegies Zwart et alll200ll l 
and the "NBODY6++" (jSpurzemlll999r i utilize analogous 
wind treatments. The w ind mass loss being cr ucial for very 
massive stellar entities ijGlebbeek et al.ll2009l ). we consider 
the effects of more detailed stellar winds separately in Sec. [4] 
Mergers among MS stars are treated in NBODY6 by 



adopting the schemes of IHurlev et~ai1 l|2002l . 120051 ). When 
two MS stars collide (a collision between two single MS stars 
or between two MS components in a binary due to eccen- 
tricity induced by close encounters and/or due to encounter 
hardening), the two main assumptions in the above scheme 
are (i) the merged object is a MS star with the stellar ma- 
terial completely mixed and (ii) no mass is lost from the 
system during the hydrodynamical process leading to the 
merger. This no mass-loss assumption is bas ed on the results 
of SPH simulations of MS-MS mergers fe.q JSills et~aT]|200ll 
who show that the mass loss is up to 5% for low-mass MS 
stars) but such studies also show that the mixing is only par- 
tial. The age of the merged MS star is assigned based on the 
amount of unburnt hydrogen fuel gained by the hydrogen- 
burning core as a result of the mixing. In case of a mass 
transfer across a MS-MS binary, the cases of the original 
accretor MS star having a radiative and a convective core 
are treated separately (included in the BSE algorithm). For 
a convective core, the core grows with the gain of mass and 
mixes with the unburnt hydrogen fuel so that the accreting 
MS star appears younger. For the case of a radiative core, the 
fraction of the hydrogen burnt in the hydrogen-burning core 
remains nearly unaffected by the gain of mass so that the 
effective age of the MS star decreases. However, in any case 
of a merger between two MS stars, an appropriate amount 
of mass is removed from the final product if the kinetic en- 
ergy of their final approach exceeds the binding energy of the 
merged MS star. Admittedly, there exists some arbitrariness 
in applying the treatments, appropriate for mergers of low- 
mass stars, to the cases involving mergers of w lOOM© stars 
which we consider in the following sections. The structure of 
these stars is very different so that it is unclear whether the 
above collisional mass-loss and mixing criteria will remain 
valid for such massive stellar mergers. Unfortunately, the 
outcomes of massive stellar collisions are currently unclear 
from theoretical studies. A relatively simplistic but very fast 
algorithm f or treating collis i ons o f massive stars (the MMAMS 
scheme) by iGaburov et~ai1 ((2008) shows < 10% mass loss 
and partial mixing in encounters of ~ IOMq stars. In the 
absence of a clear understanding, we continue to utilize the 
scheme for low mass stellar mergers for massive mergers also. 

Note that the computed Models 1-4 of Paper I consti- 
tute exactly the same cluster model of 4 different random- 
ized (discrete) realizations. Each of these models are evolved 
upto « 3 Myr which is the widely used ag e of R136 as in- 
ferred from its low mass stellar population l|Andersen et al.l 
2009). The above age estimate is perhaps the most robust 
one for the R136's stellar population. 

2.3 Additional model computations 

In the four model calculations presented in Paper I (Models 
1-4), there are typically 1-2 massive binary mergers per clus- 
ter at the very beginning of the computations that form SC 
stars. These mergers are found to happen mostly due to high 
initial eccentricities of the corresponding massive primordial 
binaries. Although these merger products are highly super- 
canonical, we consider such immediate mergers as "acciden- 
tal" since the parent binaries would have already merged 
during their pre-main-sequence evolution. However, we ad- 
here to the canonical upper limit of m max * = 150Mq for 
R136 which is a plausible assumption as the limit is deduced 
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Figure 3. Radial distances, i?sCi °f the SC stars w.r.t. the cluster's center of density that are formed in the models "3", "C5" and 
"CIO". It can be seen that while the SC stars generally form and remain close to the cluster's center, some of them form with runaway 
velocities (the monotonically outgoing trajectories; see Sec. 13 - 1 1) and they soon escape as runaway stars from the cluster. The colours of 
the lines correspond to the same stars as in Figs.[l]&[2] 



from observations of young massive star clusters. Hence, we 
consider the SC stars, that appear due to a collision be- 
tween the two members of a primordial binary at the very 
beginning of the cluster evolution (c./. Fig. [I}, as spurious 
and unrelated to the formation of the initial ZAMS stellar 
population. These are artifacts of the chosen eccentricity dis- 
tribution. The SC stars that are formed later in the course 
of the cluster evolution are, of course, considered genuine. 

Notably, the presence of 1-2 spurious (single) SC stars 
is not expected to affect the dynamical evolution of our 
model clusters in any significant way as for a massive primor- 
dially mass-segregated system like ours, there are anyway a 
large number of massive stars in the cluster's central region. 



Therefore the inclusion of a few more massive single stars 
is unlikely to provide a substantial additional effect. It is, 
of course, not impossible that two massive stars happen to 
collide to form a single SC star at the birth of a cluster, but 
it seems unphysical that it could be an eccentricity effect 
which is always the case for the above mentioned acciden- 
tal SC stars, which is why we consider them spurious. The 
conditions under which such "primordial-mergers" can hap- 
pen and their chances is beyond the scope of the present 
study. In any case, apart form their negligible contribution 
to the dynamics of their host clusters, it is unlikely to find 
any such primordially-merged star in its SC state in R136 
at its current age, due to the short SC lifetime, so that the 



Super- canonical stars in R136 7 



currently observed ones should all be dynamically formed. 
We elaborate the latter point in Sec. 14.11 

To see whether the chosen thermal eccentricity distri- 
bution in the computed models of Paper I does have a 
systematic effect on the dynamical formation of SC stars, 
we perform 5 additional computations with the same ini- 
tial configuration (but different realizations) as in Paper I 
(see Sec. 12. ip except that all the primordial binaries are 
initially taken to be circular. Of cours e, the most realis- 
tic approach would be to "eigenevolve" (|Kroupal ll995'). i.e., 
to incorporate the pre-main-sequence evolution of each of 
the primordial binaries to determine their eccentricities and 
mass-ratios at the beginning of the computations. However, 
the nature of eigenevolution is not yet known for progenitors 
of massive stars. Typically, the tidal interaction between the 
proto-massive-star members would circularize the tight bi- 
nary orbits but wide binaries would still remain eccentric. 
Currently, it is unclear in the literature what would be the 
exact outcome of eigenevolution of massive proto-stellar bi- 
naries. As we shall see in Sec. [3] if only the dynamically 
formed SC stars are considered, their population is similar 
for the calculations in Paper I and in these new computa- 
tions, i.e., the primordial binaries' eccentricity distribution 
does not play a significant role. 

We discuss the formation of SC stars from the com- 
putations of Paper I and from these new computations in 
Seed 



3 RESULTS: FORMATION OF 
SUPER-CANONICAL STARS 

Fig. [1] (top panel) shows the masses of the SC stars, Msc, 
with the cluster evolution time for the computed "Model 
3" of lBaneriee. Kroupa fc Ohl l|2012h . Each appearance of a 
line in this figure corresponds to the appearance of a SC 
star and its slow decline is due to the wind mass loss of 
the SC star as provided by the currently available stellar 
wind prescriptions in NBODY6 (see Sec. 12 . 2f) . As pointed 
out in Sec. 12.21 this stellar wind scheme is inappropriate 
for very massive stars like the super-canonical ones, as it 
substantially underestimates the mass loss and enhances the 
lifetime in the SC phase. We address this issue in Sec. [4] 

The SC stars represented by the grey lines that appear 
from the very beginning of the evolution are formed due to 
immediate mergers of highly eccentric massive primordial 
binaries which are unintended and considered spurious as 
explained in Sec. 12.31 The rest of them are of course gen- 
uine SC stars that are formed through dynamically-induced 
mergers of massive binaries (see Sec. 13. 1 P - This particular 
model formed a rather marked number of SC stars. Notably, 
a Msc ~ 240Mq SC star is formed for the first time as early 
as To = 0.7 Myr. 

Fig. Q] (bottom panel) similarly shows the formation of 
SC stars for "Model 4" of Paper I. Here, only one spurious 
SC star is formed which is only marginally more massive 
than the canonical upper limit and hence disappears soon 
due to wind mass loss. The subsequent cluster contains only 
dynamically formed SC stars. In this model also, the first 
appearance of a SC star (Msc ~ I7OM0) occurs quite early 
(To = 1.2 Myr). 

Fig.[2]shows the formation of SC stars in the additional 



models "C2", "C5" and "CIO" with all primordial binaries 
initially circular as explained in Sec. 12.31 While the total 
number of SC stars formed in these models are typically 
similar to those in the models of Paper I (c.f. Table [T]), 
none of them is formed artificially as can be seen in Fig. [2] 
The SC stars start to appear from typically To « 1 Myr 
for these models also. The remaining two of the addition- 
ally computed models produce SC stars only with masses 
marginally above 15OM0. 

3.1 Dynamical formation of SC stars 

In the above computations, we find several dynamical chan- 
nels for the formation of the SC stars. The most frequent 
way is found to be an abrupt eccentricity enhancement of a 
massive binary due to a close encounter with a single star. 
A collision and subsequent merger occurs between the com- 
ponents when their periastron separation becomes smaller 
than the sum of their radii. The single stars in the cluster's 
central region arise from dissociation of wider binaries by en- 
countering with harder binaries or by direct ionization (i.e., 
detachment of a binary due to its encounter with a single 
star of kinetic energy larger than its binding energy). This 
channel is most fruitful for the Models 1-4 as they generally 
contain more eccentric binaries. 

The next important channel is the eccentricity augmen- 
tation of the inner massive binary due to Kozai cycles with 
an outer companion. Such hierarchical triples can often form 
in the cluster center as a result of binary-binary encounters. 
Notably, a close encounter of a single star with a hard binary 
can also form a temporary tri ple system (resonance encoun- 
ters; see iHeggie &: Hun 12003). The orbital changes due to 
tidal interactions within the inner massive binary and also 
that of the outer member with the inner members and as 
well the merger process of the inner binary can often make 
the triple unstable. This causes the inner binary to become 
isolated from the outer companion in the course of its ec- 
centricity growth. In our computations, a few newly formed 
SC stars are still found bound in wide orbits with the outer 
member which soon get ionized (within a few orbital times). 
In a few cases, a SC star is also found to form when the outer 
member of a triple collided with one of the inner m embers. 

SC stars also form due to encounter hardening fHcggic 

Il975h of the orbits of appropriate binaries resulting in a 
Roche lobe overflow. The approach to the Roche lobe con- 
tact is generally assisted by the evolutionary expansion (dur- 
ing its main sequence) of the more massive bi nary member 
which therefore fills the Roche lobe first (c.f. iHurlev et al.l 
120051 ). The subsequent unstable mass transfer from the more 
massive member to the less massive one coalesces the binary. 

Some of the SC stars are found at high velocities at their 
birth. These are the ones whose progenitor binaries suffered 
substantial recoils in the close encounters that caused them 
to merge and they escape from the cluster right at their 
birth. As examples, the radial distance, Rsc, vs. t plots for 
the SC stars in Fig. [3] show such runaways in the models 
"3" and "CIO" exhibiting the monotonic outgoing trajec- 
tories. In fact, such born escapers constitute most of the 
runaway stars in the SC mass-range. This can be expected 
since once formed, it is unlikely to eject a SC star later by 
a super-elastic close encounter (lHeggiel [T97flh with a hard 
primordial binary unlike the non-SC stars for which such an 
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ejection channel is usual. The SC star, being the most mas- 
sive participant in such an encounter (the primordial binary 
members are < 150M©), is likely to get trapped in the bi- 
nary by exchanging with one of the companions and hence 
receiving less recoil. 

The rest of the SC stars initially form and remain bound 
to the cluster. Because they are among the most massive 
members of the cluster, they continue to reside close to the 
cluster's center, occasionally several of them simultaneously 
(c.f. Fig. [3J). However, a few of them can subsequently be 
ejected by dynamical recoils (only two such SC ejecta oc- 
curred from all of our computations). Notably, the ejection 
or runaway fraction of stars in the SC mass range is > 0.5 
as determined in Paper I (see Fig. 4 of Paper I). 

It is to be noted that the above dynamical channels 
are as well applicable for formation of merger products with 
m s < 150M Q . 



4 EVOLUTION OF SUPER-CANONICAL 

STARS: THE EFFECT OF STELLAR WINDS 

As already raised in Sec. 12.11 the treatment of stellar wind 
mass loss in NBODY6 is inappropriate for massive (main- 
sequence) entities like the SC stars which grossly underesti- 
mates their wind mass loss, thereby enhances their lives in 
the m 3 > 150Mq state. In Fig.[T] it can be seen that all the 
SC stars formed in "Model 3" remain super-canonical until 
the end of the computation at t ~ 3 Myr. This is as well true 
for all the substantially super-canonical merger products in 
the other computations. 

To take into account an appropriate mass-loss for the 
SC stars, we consid er the stellar evolutionar y models of ro- 
tating MS stars by iKohler fc Langerl (|2012h which encom- 
pass the SC mass range and are computed for the LMC 
metallicity. The above authors use the one-dimensio n al hy - 
drodynamic binary evolution code of iHeeer et all {2000) 
which includes direct treatments of stellar rotation, mag- 
netic field and detaile d stellar winds. The c ode has been re- 
cently improved as in lPetrovic et all (|2005l ) and lYoon et alj 
<|200fih . 

Stellar wind is the most important factor for the evolu- 
tion of the mass of an isolated ma ssive star during i t s mai n 
sequence. In the comp utations bylKohler fc Langerl l|2012l ). 
the MS stellar wind of IVink et all (|200ll ) is applied that in- 
cludes the bi-stability jump. Durin g the bi-stab i lity ju mp, a 
linear interpolation is applied as in iBrott et al.l l|201ll ). The 
mass loss rate is switched from | Vink et al.l ((2001) to that of 
iNieuwenhuiizen fc de Jager] |l99ol ) when the latter becomes 
stronger (at effective temperatures < 22000 K). When the 
surface Helium abundance becomes Y s 0.7, the star is 
considere d to enter the WR ph ase in which case the much 
en hanced lHamann et al.l (|l995l ) wind, divided by a factor of 
10 l|Yoon et al.ll2006l ). is applied. 

Fig. [4] shows the resulting mass evolution of SC stars 
with initial surface rotation velocities of Kot(0) = 50 & 
500 km s _1 . The SC stars undergo a steep mass loss as they 
become Wolf-Rayets. Because of the large mass loss in the 
WR phase, these stars become less massive than 15OM0, 
i.e., cease to be super-canonical, in rsc ~ 1-4 — 1.8 Myr 
from their zero-age. The lifetime in the SC phase therefore 
becomes much shorter with the more detailed mass loss pre- 



scription in a hydrodynamical stellar evolution calculation. 
We consider this more realistic lifetime as the lifetime of a 
SC star rather than that obtained from within NBODY6. 
We note that the SC stars formed in our computations do 
not undergo further coalescence or Roche lobe overflow with 
a secondary (although these are in principle possible) so that 
the individually obtained mass loss rates can be applied to 
them (see also Sec. [5}. It can be noted from Fig. [4] that the 
mass loss rate is only moderately affected by the stellar spin. 
The SC stars being binary merger products are indeed likely 
to form with high surface rotation velocities. 

Notably, in a dynamically active environment, a SC 
star can undergo further mass evolution due to subsequent 
merger of the SC product with other stars or its acquire- 
ment of a binary companion causing Roche lobe overflow. 
However, in our computations, we hardly find subsequent 
mergers of the SC stars or them being overflown as con- 
firmed by their continuous and smooth mass depletion (c.f. 
Figs. Q] & [2}, justifying the application of isolated stellar 
mass loss to them. 



4.1 Can the dynamically formed SC stars in R136 
be observed? 

Given that the lifetime of the SC phase (for the LMC metal- 
licity) is rsc ~ 1-4 — 1.8 Myr, how likely is it to detect the 
dynamically formed SC stars in R136 at the present day? 
Clearly, the SC stars must form at most a rsc time earlier 
to be presently visible in their SC phases. From Figs. [1] [2] 
[3] and Table [T] (see below) , it is clear that (a) SC stars are 
very likely to form at any time within 1-3 Myr cluster age 
and (b) multiple SC stars are likely to appear and remain 
bound to the cluster within a SC lifetime of rsc = 1-5 Myr, 
the latter value being taken representatively. 

These conclusions imply that in spite of the severe wind 
mass loss of the SC stars, it is still feasible that they can 
be detected at the present day in multiple numbers near the 
center of a R136-type cluster. Such SC stars are the ones 
that are formed via dynamical mergers of massive binaries 
late enough that they remain super-canonical until now. Our 
calculations indicate that the satisfaction of this condition 
is quite feasible for a R136-type cluster. For example, at the 
cluster age of 2.5 Myr, "Model 3" would retain 3 SC stars 
[c.f. Figs. [T]and |3}, C5 and C10 would retain 3 and 2 SC 
stars respectively (c.f. Figs. [5] and [3} . 

By the same argument, the above detailed evolutionary 
models of the SC stars imply that those observed in R136 
would have evolved below the canonical limit by the present 
day had they been primordial, i.e., had they been formed 
at the cluster's birth, thereby violating the canonical up- 
per limit, assuming th e widely used « 3 Myr age of R136 
If Andersen et al\\200!\) . This strongly supports our dynam- 
ical formation scenario of the observed SC stars in R136 
which causes them to form sufficiently later so as to make 
them visible during their SC phases at the present day. This 
is well supported by our computations. 

Admittedly, there are several drawbacks and incom- 
pleteness in the above line of arguments. First, t he de tailed 
stellar evolutionary models of Ko hler fc Lan gcr (2012|) have 
been considered independently of the dynamics of the model 
clusters. The newer stellar evolutionary models would result 
in stronger wind mass loss from the cluster core and would 
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Figure 4. Evolution of the mass of initially SC stars as computed bv lKohler &: Langerl J2012f ) for initial surface rotation velocities of 
Vrot (0 ) = 50 & 500 km s — 1 and metallicity appropriate for the LMC. Followed by an initial moderate mass loss, as given by the lVink et al.l 
(2001) stellar wind, a SC star undergoes a substantial mass loss as it enters the WR phase. The black horizontal line highlights the 
canonical upper limit of 15OA-/0. It can be seen that the stellar mass remains super-canonical until rs 1.4 — 1.8 Myr from the ZAMS. See 
text for details. 



contribute to the dilution of the cluster's central potential to 
a larger extent than t hat from the NBQDY6's native stellar 
evolutionary scheme (|Hurlev et al ] |200d . l2002h . However, in 
our computations the cluster's core, where the binary merg- 
ers occur, is populated by a large number of hard massive 
binaries due to our initial conditions (Sec. I2.1|l . which are 
all in accordance with observations of young star clusters 
(see Sec. 12.11 and references therein). In that case, the dy- 
namics of the cluster's central region will be dominated by 
the energy generated in t he super-elastic binary-single and 
binary-binary encounters (Hcggic 1975; Hcg gie fc HutJ[2003l ) 
and the stellar mass loss would not affect the dynamical en- 
counter rates in the cluster's center significantly. 

It is true that the mass loss from the individual com- 
panions of a binary would result in expansion of the binary 
orbit and hence dilute i t s hard ness. However, as shown by 
iFuiii fc Portegies ZwartJ (|201ll), this effect is counteracted 
by the dynamical hardening (|Heggieill975l ) of the binaries 
so that t hey remain hard through o ut. In the N-body calcu- 
lations of IFuiii fc Portegies ZwartJ (|201ll ). even just a single 
binary survives to remain hard in their model cluster's core 
and dominates the energy generation in the cluster's core 
(until its ejection). With a large number of binaries, as in 
our case, the encounter hardening and the corresponding 
energy generation would be substantially stronger as the 
binaries provide much larger encounter cross sections than 
sing le stars due to t heir much bigger geometrical size (see 
also IFuiii et al.112012]) . Noticeably, even after including stel- 



lar wind mass loss, the comput ed single-star-only models of 
IFuiii fc Portegies ZwartJ (|201ll ) undergo deep core-collapse 
within ~ 3 Myr, i.e., the effects of dy namical relaxation 
and negative specific heat (|Spitzenll987r ) themselves domi- 
nate over the wind mass loss. Our models, of course, do not 
core collapse due to the substantial energy release from the 
central primordial binaries. 

Another drawback of the above arguments is that they 
are based on stellar evolution models that begin from ZAMS 
with normal composition. The SC stars are merger prod- 
ucts of O-stars that are evolved from their zero age and 
therefore would be He-enriched and substantially rotating. 
While we do consider rotating models, the initial He frac- 
tion is take n to be that appropri ate for the LMC, viz., 
Y = 0.2562 (|Kohler fc Langeijl2012T l. To obtain a basic esti- 
mate of the He-richness of the merged product, let us con- 
sider two WOMq stars that merge to form a typical 200Mq 
SC star at a typical cluster age of 1.5 Myr. Each of the ini- 
tial stars (before merger) would possess a H-burning core of 
M c « 8OM0 and at the above merging age their core He 
fraction Y c ~ 0.5 (i.e., « 40Mq He per star) according to 
the above newer models (Karen Kohler; private communi- 
cation). Assuming an ideal complete mixing and no mass 
loss during the merger process, the new He abundance of 
the merged product would be Y ~ 0.4 which would also be 
the rejuvenated surface value Y s ; for a partial mixing the 
surface abundance would be lower. While this Y s is large 
enough to drive a stronger wind from the newly formed SC 
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star compared to that for its ZAMS value Y 3 m 0.25, it is 
still moderately above the ZAMS value and less than the 
abundance at which the WR winds turns on (Y a ^ 0.7; see 
Sec. [4]). Therefore, a drastically large wind from the newly 
formed SC star is not expected. 

The above discussions imply that none of the short- 
comings or uncertainties of our approach dictate to the im- 
possibility of observations of SC stars in R136; rather such 
observations are still feasible. 



5 DISCUSSION AND SUMMARY 

Table [1] summarizes the SC star formation in the computed 
models of Paper I and in the new computations. Our general 
conclusions from the computations are as follows: 

• Formation of SC stars due to dynamically induced 
mergers of massive binaries are common in a R136-type 
cluster. In most of the computations we find multiple SC 
stars formed within < 3 Myr which is appropriate for R136 
I Andersen et al.ll2009l ; see Figs. CO & 1J. Most SC stars are 
formed single while a few of them initially form in wide bi- 
naries. 

• The SC stars typically begin to appear from To ^ 1 
Myr cluster age or even earlier (c.f. Table [l}. They tend to 
form with equal likeliness over a cluster age of 1 - 3 Myr. 

• The most massive SC star formed in a given computed 
model is typically close to M max « 2OOM0 and the most 
massive one formed is « 25OM ("Model 3"; c.f. Table CO). 
The cluster age T max corresponding to the formation of the 
most massive SC star is typically well within < 3 Myr (c.f. 
Table [TJ . 

• Multiple SC stars are occasionally found to exist simul- 
taneously near the cluster's center which are bound to the 
cluster, over a representative SC phase lifetime of tsc = 1.5 
Myr, within < 3 Myr (c.f. Table CQ) cluster age. This, along 
with the second point above, implies that it is quite plau- 
sible that R136 harbours multiple SC stars at the present 
day. 

• Some SC stars are formed with runaway velocities and 
escape (c.f. Sec. 13.11 Fig. [2]). 

These conclusion s conform with the observations by 
ICrowther etHI l |201Ch who found 4 SC stars in R136 in the 
initial mass range of 165 — 32OA/0. In our computations, 
SC stars of up to w 250M© are formed which is consistent 
with the observations taking into account the large uncer- 
tainties in the stellar evolution models in this mass range. 
Also, only one of the models ("Model 3") has up to 4 SC 
stars simultaneously bound to the cluster (over a 1.5 Myr 
period) although other models also host multiple SC stars 
close to the cluster's center. Some models however contain 
only one bound SC star within 3 Myr (c.f. Table CQl. 

At this point, it is worthwhile to note that our above 
conclusions depend somewhat on the age of the R136 clus- 
ter which is, of course, not fully settled yet. Whi le the bulk 
of R136 is « 3 Myr old (|Andersen et al.l 120091') the high- 
mass stellar population mig ht be younger ( Massev fc Hunter! 
1 19981 ; Ide Koter et al.l fl998j ) . In our computations, the SC 
stars typically begin to appear from ~ 1 Myr cluster age 
and several of them appear by « 2.5 Myr. This time-frame 
is quite consistent with the possible < 2 Myr age of R136's 



massive stellar population given the wide uncertainties in 
the age of such a young cluster particularly for the mas- 
sive end of the IMF. There has so far been no direct esti- 
mate of the age of the massive end of the stellar population 
of R136 and the above age-limit is based only on compar- 
isons with the stellar po pulations of Car OBI and NGC 
3603 (jCrowther et alj|2010h . Given that the stellar IMF o f 
R136 continues to maintain the canonical law (|Kroupall200ll ) 
from the low to the h i gh mass range (I.IM q — I2OM0; 
iMassev fc Hunted Il99gj ; lAndersen et alj I2009TI . we find it 
more natural to consider that the whole stellar population 
of the R136 cluster has formed in a single starburst event 
whithout a significant age spread. The issue can, of course, 
be more resolved with better understanding of the evolu- 
tion of ve ry massive s t ars an d their winds. Finally, a recent 
study by IChini et all |2012f ) shows that the O-star binary 
distribution can, in fact, be even harder (i.e., more bound 
or tighter) than that considered in our computations (see 
Sec. 12. ip . This would lead to the formation of SC stars even 
earlier which constitutes an important future study. 

It is important to remember that the formation of 
SC stars at early ages in the above computations is fa- 
cilitated by the adopted initial complete mass segrega- 
tion. This condition subjects the massive binaries to strong 
dynamical encounters from the beginning of the cluster 
evolution. Note however that primordial mass segregation 
is inferred to have been true for several Galactic globu- 
lar clusters jBaumgardt et alJl2008l : lMarks fc Kroupa 

b 



Strader ct al. 2011) and open clusters (ILittle fair et al 



2010; 



2003; 



Chen et al.1 120071 ). Notably, for a completely unsegregated 
model, it would take 10 Myr for the massive binaries to 
segregate to the cluster's center which makes the early for- 
mation of SC stars unlikely. However, the mass segregation 
timescale shortens substantially with increasing compact- 
ness of the cluster. In this context, an important outlook 
would be to study the formation of SC stars with varying 
initial compactness and degree of primordial mass segrega- 
tion. 

The drawbacks of the computed R136 models, as dis- 
cussed in Paper I (see Sec. 4 of Paper I), naturally carry 
over to the present analyses. These limitations however do 
not crucially affect SC formation. In particular, the trun- 
cation of the binary distribution at m 3 w 5Mq does not 
influence the mergers of the most massive binaries, the lat- 
ter being much more centrally concentrated due to mass 
segregation. Also, the exact period distribution of the O- 
star binaries is not instrumental for the formation of the 
SC sta rs as long as O-star s are largely found in tight bi- 
naries (|Sana fc E vans 2011). Notably, our adopted range of 
the orbital period for the O-star binaries is similar to that 
reported by I Sana fc Evans! (|201lh (see Paper I). The forma- 
tion of SC stars is generally similar for models with initially 
thermally distributed eccentricities (Paper I) and for the 
newer ones with initially circular binaries (except for the 
spurious SC stars in the former models). This implies that 
the initial eccentricity distribution of the massive binaries 
does not crucially influence the formation of SC stars. 

The most consistent and accurate way to address the 
above discussed concerns regarding the effects of a stronger 
wind mass loss due to the newer stellar evolutionary models 
is to incorporate the latter in a direct N-body code which 
is well beyond the scope of this paper. No such direct N- 
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Table 1. Table summarizing the formation of SC stars in our computed models. The de- 
scriptions of the columns are as follows: Col. (1): model ID, Col. (2): time To at which a 
SC star first appeared in the model, Col. (3): birth mass Mo of the first-comer SC star, 
Col. (4): birth mass M max of the most massive SC star formed in the model, Col. (5): for- 
mation time, T max , of the most massive SC star, Col. (6): maximum number of SC stars, 
A/" max ,in, that remained simultaneously bound to the cluster over a rsc = 1.5 Myr period 
within < 3 Myr cluster age, Col. (7): total number of SC stars, A/tot > formed over the 
computation (all of them do n ot necessarily appear simultane ously or remain bound to the 
clus ter). For Models 1-4 (from lBaneriee. Kroupa &: Ohll2012l ) the "spurious" SC stars (see 
Sec. 12.31) are excluded and only those SC stars that are formed later in these computations 
are considered. 



Model ID 


T (Myr) 


M (M ) 


Mmax (M ) 


T max (Myr) 


■^max,in 


Mot 


1 


2.6 


193.9 


193.9 


2.6 


l 


2 


2 


2.0 (3.0)° 


155.2 (181.4) Q 


181.4 


3.0 


i 


2 


3 


0.7 


236.8 


246.0 


1.5 


4 


5 


4 


1.2 


172.5 


206.2 


2.6 


1 


2 


C2 


1.4 


220.6 


220.6 


1.4 


1 


2 


C5 


1.3 


224.0 


224.0 


1.3 


3 


3 


C10 b 


1.2 (2.1) a 


152.4 (162. 5) a 


225.9 


2.2 


2 


4 



a The first-comer SC star's mass is too close to 150Mq. The time and the mass corresponding to the next 
SC appearance is also shown in the parentheses. 

^ The remaining two of the additionally computed models (with initially only circular binaries) produce 
SC stars only with masses marginally above 150Mq. 



body code exists at the present time. Furthermore, there are 
substantial uncertainties in the physics of mergers of mas- 
sive stars at the present time. Therefore, the conclusions 
in this work are the best one can draw given the current 
technical limitations. Perhaps it would be possible to do 
such work in future with a distributively computing, highly 
modular N-body calculation framework such as "MUSE" 
l|Portegies Zwart et alj|2009h . In spite of all these uncertain- 
ties it is worthwhile to note that if one focuses to the par- 
ticular case of R136, then (a) the existence of SC stars is 
strongly supported by observations and (b) considering the 
widely used « 3 Myr age of R136, the SC stars must have 
formed later than the birth of R136. Our computations are 
consistent with the formation of SC stars through dynam- 
ically induced mergers of massive binaries. While there are 
technologically-limited drawbacks in our present analysis, 
as elaborated in Sec. 14.11 none of these shortcomings dic- 
tate to the non-observation of SC stars in R136. Therefore, 
we can say that we have justified the hypothesis of Sec. [1] 
through our detailed N-body modelling of R136 with initial 
conditions chosen in accordance with observations of young 
star clusters. In other words, from our realistically modelled 
computations of R136-like star clusters, and from our anal- 
yses as presented above, it can now be more definitely con- 
cluded that the observed super-saturated stellar population 
in R136 does not imply a violation of the canonical stellar 
upper mass limit near m max * = 150Mq. 
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